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Abstract. We propose three implementations of the Gauss sum factorization 
schemes discussed in part I of this series [S Wolk et al., preceding article]: («) a 
two-photon transition in a multi-level ladder system induced by a chirped laser pulse, 
(ii) a chirped one-photon transition in a two-level atom with a periodically modulated 
excited state, and (iii) a linearly chirped one-photon transition driven by a sequence 
of ultrashort pulses. For each of these quantum systems we show that the excitation 
probability amplitude is given by an appropriate Gauss sum. We provide rules how to 
encode the number N to be factored in our system and how to identify the factors of 
N in the fluorescence signal of the excited state. 
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1. Introduction 



In [T] we have shown that Gauss sums are excellent tools to factor numbers. Throughout 
this article we have concentrated on the general principles and the underlying 
mathematical foundations. However, we have not addressed physical implementations of 
Gauss sums. The present article complements |T] and proposes three distinct realizations 
of Gauss sum factorization. 

Several experiments pHTT] have already successfully demonstrated factorization 
with the help of Gauss sums. However, all of them have implemented the truncated 
Gauss sum 

m=0 ^ ' 
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but none of the Gauss sums introduced in pQ . Moreover, in all these experiments, except 
the recent one using the Michelson interferometer, the ratio N/£ had to be precalculated. 
In contrast the three implementations of Gauss sum factorization proposed in the present 
paper, encode the number to be factored and the trial factor i in two independent 
variables. As a consequence, the ratio N/i does not to be precalculated. 

The first suggestion utilizes a two-photon transition in an equidistant ladder system 
driven by a chirped laser pulse. Since in this system the excitation probability amplitude 
is given by the continuous Gauss sum 



discussed in [12j, it can factor numbers. 

In addition, we suggest two alternative approaches based on one-photon transitions 
in a laser-driven two-level system. Both rely on quantum interference of multiple optical 
excitation paths, but lead to different quadratic phase factors. 

We consider a two-level system with a permanent dipole moment in the excited 
state. A cw-microwave field interacting with this dipole modulates the energy of the 
excited state and induces in this way an equidistant set of sidebands. The probability 
amplitude for a one-photon transition caused by a chirped laser pulse is the sum over 
all possible excitation channels involving one optical photon as well as multiple quanta 
of the microwave field. The Gaussian nature of this sum arises from the quadratic chirp 
of the laser pulse. 

The second approach uses a multi-pulse excitation of a two-level system with a 
linearly chirped resonance frequency. The probability amplitude of excitation is a sum 
of contributions arising from each pulse. The quadratic phase dependence characteristic 
of a Gauss sum originates from the linear variation of the resonance frequency of the 
two-level system. 

Throughout the article we neglect spontaneous emission since the interaction time 
with the pulses is much shorter than the decay time of the atomic level. This fact allows 
us to describe the system by the Schrodinger equation rather than a density matrix. 

Our article is organized as follows: Since we rely heavily on the physics of chirped 
laser pulses we first summarize in Sec. |2] the basic elements of this branch of optics. In 
Sec. [3] we then recall the main results of Ref. [12] and, in particular, the expression for 
the excitation probability amplitude corresponding to a chirped two-photon transition in 
an harmonic ladder system. Starting from this formula we demonstrate that this system 
allows us to factor numbers, and that in principle a single realization of a factorization 
experiment can be employed to reveal the factors of another number. 

Section m provides the basic elements of Sees. [5] and [6] by engineering a one-photon 
transition. As a second physical system implementing Gauss sums we investigate in 
Seels] a two-level system with a permanent dipole moment driven by a microwave field. 
A chirped laser pulse interacts with the so-engineered Floquet ladder. The resulting 
probability amplitude is again given by a Gauss sum. Section |6] is devoted to the 




m=-M 
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analysis of the third system which features a hnear variation of the resonance condition 
of a two-level atom exposed to a pulse train. This system provides an alternative way to 
factor numbers based on Gauss sums. After a comparison of both factorization schemes 
in Sec. [7] we conclude in Sec. [8] by a summary of our results and a brief outlook. 



2. Chirped pulses: Essentials 

Throughout the article we take advantage of the technology of chirped laser pulses. For 
this reason we briefly summarize in the present section the key ideas and formulas of 
this field. 

A chirped laser pulse 

Eit)^So[e~'-''fit) + c.c.] (3) 

consists of an amplitude Sq, a carrier frequency ul and a pulse shape function 

1 



fit) = fo exp 



-^(Ac./o)V 



(4) 



Here we have introduced the complex-valued amplitude 



and the dimensionless parameter 

a = Aw V" (6) 

represents the second order dispersion. Moreover, Au denotes the bandwidth of the 
pulse and 0" = d^0(c<;)/dco'^ is a measure of the quadratic frequency dependence of the 
phase of the laser pulse. 

When we substitute ([5]) into the exponent of the Gaussian in (jl]) we find that the 
pulse shape 



of a chirped pulse consists of the product of a real-valued Gaussian and a phase factor 
whose phase is quadratic in time. 

Since the instantaneous frequency 

d / aAcj^ ,\ aAw^ , , 

iy(t) = — — = -t (8) 

^ ^ dt V2(l + a2) J 1 + ^ ' 

of the pulse is the derivative of this phase with respect to time, the frequency changes 
linearly in time as the pulse switches on and off. For a positive value of 0" we find an 
increasing frequency, whereas a negative (j)" corresponds to a decreasing frequency. 
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Figure 1. Model of ladder system. The ground state | g) is connected by a two-photon 
transition to the excited state |e). We include an harmonic manifold oi D = M' + M+1 
intermediate states \m) with M' <m<M which are shifted by the offset 6m = A 
with respect to the central frequency wo- The offset of the central state with to = is 
i5, whereas neighboring states in the harmonic manifold are separated by A. 



3. Chirping a two-photon transition 

In Ref. [12] we have considered a two-photon transition in the ladder system of Figured] 
driven by a chirped laser pulse. In the weak field limit, the probability amplitude to 
be in the excited state results from the interference of multiple quantum paths each 
contributing a quadratic phase factor. Since the population in the excited state has 
the form of a continuous Gauss sum we can use this observable to factor numbers as 
suggested [1] in part I of this series. 

We first briefiy summarize the essential results of Ref. [I2j. Then we turn to a 
demonstration of the factorization capability and bring out the physical origin of the 
scaling property derived in Ref. [1] from mathematical arguments. 



3.1. Brief review of the model 

We consider a quantum system with a ground state \g) and an excited state |e) separated 
by an energy 2fvjjQ. In the neighborhood of the midpoint of the energy difference we 
assume to have a manifold of equidistant energy levels as shown in Figure [1] Their 
off-set 

5m = 5 + mA (9) 

with respect to the central frequency Uq is the sum of the off-set 6 to the central state 
and integer multiples m of the separation A of two neighboring states of the manifold. 

At this point it is important to note that this decomposition of 6m is not unique. 
Indeed, we could have chosen a " central" level which is different from the one indicated 
in Figure [1] This choice would have changed the integers m. This ambiguity in the 
labeling of the states is the deeper physical origin of the scaling property of the Gauss 
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sum already mentioned in Ref. [T]. 

Next we drive this ladder system by a chirped laser pulse of the form given by (|3]). 
In second order perturbation theory the probability amplitude 

^(TPT) ^ ^ng^ (IQ) 

to be in the excited state after such a two-photon transition is given [12], apart from 
the phase factor exp(z7) by the Gauss sum 

M 



2vri(m + — K 



^ Wm exp 

m=-M' 

discussed in part I of this series. Here the variable 



m 



TT 



'111 



(12) 



is expressed in terms of the parameters 6 and A of the harmonic manifold of intermediate 
states and the dimensionless chirp (p". Hence, by varying the chirp we can tune C,- For 
this reason we call the variable ^ the dimensionless chirp. 
The number 

25 



N = 



(13) 



to be factored is represented by the ratio of the two characteristic frequencies of the 
ladder. 

The weight factors 

2" 



erfc ( i-T— Vl — ia ] exp 



5m 



(14) 



contain the complementary error function 

oo 

2 



erfc(z) = 



du e"" 



of complex argument z, and the abbreviation 



(15) 



(16) 



2 Aoo^ 

involves the Rabi frequencies flmg and flem connecting the ground and the excited state 
with the intermediate states, respectively. 



3.2. Factorization 

We are now in the position to discuss our factorization scheme. For this purpose we 
first address the experimental requirements and limitations and then present an example 
demonstrating the capability of this system to factor numbers. 
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3.2.1. Requirements The probability amplitude c^e^'^\ ( fTOjl . of populating the excited 
state is of the form of a continuous Gauss sum discussed in part I of this series. Thus the 
dependence of the population the excited state on the dimensionless chirp ^ 

can be employed to reveal the factors of an integer N . For this purpose we encode 
according to f|T3|) in the parameters 5 and A of the harmonic manifold of intermediate 
states. In order to apply the factorization scheme to a broad range of numbers iV, we 
require control over 5 and A. 

We emphasize that the equidistant spacing within the harmonic manifold is essential 
for obtaining the Gauss sum and for our factorization scheme. Moreover, the dimension 
D = M' + M + 1 of the intermediate levels has to be adapted to the number to 
be factored. The larger N the more intermediate states are required for a meaningful 
signal. For the factorization of we require N < D < 2N where the lowest quantum 
number is bound by M' > N/2. 

We can read out the population of the excited state by its fluorescence. In Ref. P 
we have formulated a rule for determining the factors of N: If the signal shows a distinct 
maximum around the integer ^ = £ then £ is a factor of A^. In order to resolve the signal 
in the vicinity of candidate prime numbers, we require sufficient stability of and accuracy 
in the chirp (p". 

According to Ref. [1] we also need to impose a restriction on the weight factor Wm, 
f fT4|) of the contribution arising from the quantum path through the mth intermediate 
state: Indeed, Wm must be slowly varying as a function of m in order to ensure 
that no specific excitation path is favored or discriminated. As a result the dipole 
matrix elements associated with all possible transitions should be of the same order of 
magnitude. 



3.2.2. Example Next we present numerical results for an artificial ladder system 
consisting of D intermediate states. Here, we have made the idealized assumption that 
the Rabi frequencies associated with the sequential path l^f) — )• |m) — )■ |e) are identical, 

that is Qern^mg = COUSt.. 

In Figure [2] we display the population |ce"^^"^'*p of the excited state for the number 
A^ = 15 = 3 ■ 5 as a function of the dimensionless chirp ^. For a single intermediate 
state, there exist |T2] several interfering quantum paths only if ^ < 0. However, 
for an equidistant manifold, the population is symmetric with respect to ^, that is 

The insets magnify the signal in the vicinity of trial factors = 2, 3, 5 and 7. 
Following our criterion of a dominant maximum at an integer indicating a factor we 
clearly identify from the distinct maxima at ^ = 3 and 5 the factors of A^ = 15. In 
contrast, the signal at non-factors such as ^ = 2 or 7 does not show any characteristic 
features. 
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Figure 2. Factorization ofA^ = 15 = 3- 5 with the help of the population Ici"^'^"^' p of 



the excited state given by ((TU)) and ([TT]) after a chirped two-photon transition through 
the intermediate levels of the ladder system of Figure [T] In the center we provide an 
overview over the complete signal as a function of the dimensionless chirp ^. The 
insets magnify the signal in the vicinity of candidate prime factors. Pronounced 
maxima at the prime factors ^ = 3 and 5 are clearly visible. In contrast, at non- 
factors ^ = 2 and 7 the signal does not exhibit any peculiarities. Parameters are 
S = 0.0225fs"\ A = 0.003fs"\ Aw = 0.1525fs"\ a = -10824 and 2M + 1 = 23. 



3.3. Factorization by rescaling 

In part I of this series |T] we have aheady shown that the reahzation of the Gauss sum 
Sn for is sufficient to reveal the factors of another number A^'. Whereas the proof 
presented in Ref. [1] had rehed on mathematical arguments we now use the freedom in 
labeling the levels of the equidistant ladder system to verify this surprising feature. 

For this purpose we recall from f|T3l) that N is encoded in the ratio 25/ A described 
by the offset 5 of the reference level |m = 0) and the level spacing A. However, when 
we choose a different reference level \k) with the associated offset 

5' = 5 + kA, (17) 

the number to be factored is 
26' 

N' = — = N + 2k. (18) 
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With the help of the definition (1131) for the number to be factored in terms of the 
offset 6 and the detuning A we find from f|T2l) that the dimensionless chirp 

e = 0" (19) 

71 

is proportional to N. 

As a result, the dimensionless chirp 

e' = N'—(t)" (20) 
vr 

corresponding to A^' is related to ^ and by the scaling transformation 

r = (21) 

(TPT) 

Therefore, we can factor the number A^ by analyzing the signal cl , which was 
recorded in its dependence on C, for A^, using the new scale ^' = (N'/N) ■ C,- 



4. Engineering a one-photon transition 

In the preceding section we have used a chirped two-photon transition going through an 
equidistant ladder system to factor numbers. Unfortunately, the requirements on the 
system are rather stringent and it is hard to identify quantum systems with such an 
arrangement of levels. For these reasons we now study more elementary models based 
on a two-level system with a ground state \g) and an excited state |e) separated by an 
energy hwo. 

We assume the excited state to have a permanent dipole moment pee which interacts 
with a time-dependent modulating field = Ejji{t). In the following sections we 
consider two cases: (z) a sinusoidal time dependence manifesting itself in a periodic 
modulation of the excited state, and (ii) a quadratic chirp reflecting itself in a linear 
shift. 

In addition to Em we have a time- dependent weak driving field E^ = E^{t) causing 
transitions between the ground and the excited state. Depending on the two cases E^ 
is either a single chirped laser pulse, or a sequence of pulses. 

This arrangement corresponds to the interaction Hamiltonian 

V = -PeeEm{t)\e){e\ - {pge\e){g\ + c.c.) E^{t) (22) 

where pge denotes the dipole moment of the two-level transition. Here we have assumed 
that the frequencies of i?m and E^ are clearly separated. Hence, E^ only acts on the 
excited state and E^ only on the transition. 

In the interaction picture the equations of motion for the probability amplitudes 
Cg = Ce(t) and Cg = Cg(t) to be in the excited and ground state read [H] 

i^ Ce(t) = -n^eit) Ce(t) - ^^3e(^)e'"«* Cg{t) (23) 

or 

i^c,(t) = -fie,(t)e-^""*ce(t). (24) 
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Here we have defined the time- dependent Rabi-frequencies 

^ee{t) = Pee Emit) / h and Elegit) = peg E^{t)/h (25) 

associated with the two electric fields Em and E^, respectively. 

In order to solve (123!) and (12^ we first recall that the strong field £"111 which causes 
the modulation of the excited state appears through fl^e in (1221) and multiplies Cg. Only 
the weak driving field E^ which enters (|22l) via Qge and multiplies Cg induces transitions. 
For this reason it suffices to describe this process by perturbation theory of first order. 

At time to the two-level system occupies the ground level, that is Cg(to) = 
1 and Ce(to) = 0. In the weak field limit the probability amplitude for the excited 
state does not change significantly under the action of a weak chirped pulse which yields 
Cg{t) ^ 1. Hence, (l23l) reduces to the inhomogeneous differential equation 

'Jt ""'^^^ ~ " fi,e(t)e^"°* (26) 

where the interaction with the chirped laser pulse acts as an inhomogeneity. 
It is easy to verify that the solution of (!26|) reads 

t 

Ce{t) = ie^^(*) J dt' exp[-i/3(t')]exp(iwot') ^geit') (27) 
to 

where we have introduced the phase 



t 



(3{t) = I dt' fiee(t') = ^ I dt' E^t'). (28) 



h 

to to 

When we substitute the electric field 

E^{t)=S4e-'''^'h{t) + c.c.] (29) 

of amplitude Sd, carrier frequency ojl and envelope h = h(t) into the solution (E7|) for 
the probability amplitude Ce we find in rotating wave approximation 

t 

Ce{t) = ingee'^^''^ I dt' exp[-i/3(t')] exp(i(5t') h{t'). (30) 



to 

Here we have defined the time-independent Rabi frequency 

VLge = pge£d/h (31) 

associated with the electric field £d of the transfer pulse and the detuning 5 = Uq — ujl 
between the atomic and the carrier frequency. 
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Figure 3. Engineering the Floquet ladder. We consider a two-level system with 
the ground state \g) and the excited state |e) separated by the energy huJo. The 
excited state is modulated by a strong sinusoidal field Em = Em{t), (15^ giving rise 
to equidistant sidebands separated by HA. The one-photon transition is driven by 
a chirped laser pulse — E^{t), p9|. characterized by a linear variation of the 
instantaneous frequency to = uj{t). 



5. Floquet ladder 

So far we have neither specified the modulating nor the driving field. In the present 
section we consider a sinusoidal modulation of the excited state by a strong cw field 
creating a set of equidistant sidebands, very much in the spirit of the harmonic manifold 
of Figure [TJ Moreover, we include a weak chirped laser pulse driving the transition, 
that is the envelope h = h{t) of E^i is given by f{t) of dl]). Figure [3] summarizes this 
engineering of the Floquet ladder. 



5.1. Excitation probability in the weak field limit 

We now evaluate the probability amplitude ci^^^ given by (!28|) and (130|1 for the excitation 
of the Floquet ladder. Here we proceed in two steps: we first include the modulation 
and then calculate the remaining integral for the case of a chirped pulse. 



5.1.1. Sinusoidal modulation In the case of the modulation field 

Emit) = J'o cos(At + y?), (32) 

with period 27r/A, amplitude J-q and phase ip, the time-dependent phase /3 = [3{t) 
defined by (!28|) takes the form 

= Ksin( At + (33) 

Here we have chosen the lower integration limit to = (^o^r — (/?)/ A and have introduced 
the dimensionless ratio 
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of the time-independent Rabi frequency = pee -^o/^ ^^nd A. 
When we apply the generating function [T3] 



exp (i K sm ( 



K] e 



(35) 



of the Bessel function J„ to evaluate the phase factor in the integrand of (!30l) we find 
the probability amplitude 

cf ^)(t) = ifi,ee^^« Yl '^"('^) (36) 

n 

to be in the excited state. Here we have interchanged the order of integration and 
summation and have introduced the integral 

t 

(37) 



hn{t) = j dt' exp (i 5„t') h{t') 



to 



with the offset 



5r, = 5 — nA 



(38) 



of the n-th satellite of the excited state in the splitted manifold. 

Hence, the modulation of the excited state causes equidistant sidebands and Jri(/t) 
determines the weight of the n-th sideband. In the language of Floquet theory we have 
replaced the time dependent Hamiltonian by an infinite dimensional Floquet matrix. 



5.1.2. Chirped pulse So far our calculation is valid for an arbitrary pulse shape 
h = h{t). We now perform the integration for the case of a chirped pulse where the 
envelope h is given by the complex- valued Gaussian / defined by (jl]). 

Since we are interested in times after the pulse has interacted, that is for 
VT+o^/Aw ^ t we extend the upper and lower limits of the integration in (137|) to 
+00 and — oo, respectively, which yields 



hnit) = hn = h dt e" 



-|(Aa;/ot)2+i(5„t 



or 



'111 



exp 



(39) 



(40) 



When we recall the definition (|5]) of /o we can decompose this Gaussian into a real- valued 
one and into a quadratic phase factor, that is 

2" 



^271 



Alo 



exp 



5n 

Au 







exp 





(41) 



Here we have also made use of (El). 
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In the expression for hn given by (14T1) the off-set (5„ which depends hnearly on n enters 
quadratically. Hence, hn contains quadratic phase factors. In the present section we 
cast the probability amplitude c^f^"^ given by ( 136|) into the form of a Gauss sum which 
allows us to factor numbers. 

With the help of ( BTj) we find from fl36l) the formula 



27r 



JFL) _ -r) pi/3(t) 



exp 



-imp 



exp 



Next we recall the definition (l38l) of 5ri and express it by 
5n _ n- 6/ A _ n- N 
Aw ~ Aw /A ~ An ' 
where in the last step we have introduced the abbreviations 
6 , . Aw 



N 



A 



and An 



A 

Likewise, we obtain from f l55]) the identity 
Sl = 5^-{n-n^^^)25A. 



When we substitute 
SA ... 



and fH^ into 



define the dimensionless chirp 



TT 



and recall the definition (jH]) of we arrive at the probability amplitude 



:f^)(t)=Ar(t)5^«;„exp 



ivr n — 



2n) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



to be in the excited state in the Flo quet- ladder scheme. Here, we have used the 
abbreviation 

.a 



Ar(t) = 27ri^e'^We^'^'^"/2 



together with the weight factor 

1 



^An 



exp 



1 /n~N 

2 V An 



JjK)e~ 



■imp 



(48) 



(49) 



In order to reduce the influence of the Bessel function J„ in Wn we adjust the 
modulation index 

determined by the microwave field, fl32l) . such that iy„ is slowly varying as a function of 
n. For this purpose we recall [I3] the asymptotic expansion 

y 



TT TT 

COS ( z — n 

vr^ V 2 4 



(51) 
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of the Bessel function in the hmit of large arguments n <^ z. Indeed, we find for the 
choice 



71 

K = ZTTS + — , 
4 

where s is a large integer the approximation 



JJk) = \ I — 

7YK 




l)'" for n = 2m 
for n = 2m + 1 



(52) 



(53) 



~ ( F f ) 

Thus all weight factors Wn with odd index n vanish and the probability amplitude Ce 
given by ( l47Il reduces to 

.2 



Wm exp 



m 



-27ri(m - — )e 



where 



1 1 
Wm. = — : 7= exp 



-2 



m-N/2\ 



An 



(54) 



(55) 



We note that for the choice of = 7c/2 the phase factor in Wm is unity and the 
only m-dependence left results from the Gaussian. For an appropriate choice of An, 
which according to (jH]) is determined by the bandwidth Au of the chirped pulse, this 
Gaussian is slowly varying. 

When we compare the probability amplitude Ce to be in the excited state given 
by (!54|) to the generic representation 

.2- 



S{^; A,B) = '^ Wm exp 



, m m, ^ 



(56) 



of a Gauss sum discussed in part I of this series [T] we find 

cf^)=Ar5(e;-l,iV). (57) 

Since the fluorescence signal of the excited state is proportional to the population |ce^^'' p 
in this state, it is proportional to the Gauss sum \S{C,', —1,N)\'^. 



5.3. Factorization 

In order to gain information on the factors of an appropriately encoded number N 
we present now two schemes: The first one requires a continuous measurement of the 
fluorescence signal as a function of the dimensionless chirp ^. For the second one it 
suffices to acquire the fluorescence signal at integer values ^ = L 
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5.3.1. Continuous tuning of chirp According to Ref. £ is a factor, or a multiple of 
a factor, of if — 1, A^) given by f l56|) shows a pronounced maxima aX ^ = I. In 
Figure m we present numerical results for the factorization ofA^ = 21 = 3- 7 employing a 
truncated Floquet ladder covering 2M + 1 = 78 harmonics for two choices of the relative 
phase 99. Here, we display |ce'^^'*p based on ( HTll . Again we indicate candidate prime 
factors by vertical lines. 

For if = tt/2 we clearly identify from the insets on the top the prime factors i = 3 
and 7. In contrast, the signal does not show any peculiarities at non-factors such as 
i = 2 and 5 as shown by the insets at the bottom. 

For = the phase factor e'™'^ = (—1)™ leads to oscillatory weight factors and 
our criterion of finding pronounced maxima at factors of A^ is not applicable here. 
Nevertheless, the signal |ce'^^^p still contains information on the factors of A^. Indeed, 
the corresponding signal shown at the bottom of Figure H] vanishes at the factors 
i = 3 and 7 as indicated by the insets on the top but displays no pecuharities at non- 
factors such as i = 2 and i = 5 depicted at the bottom. 

5.3.2. Discrete values of chirp Next we present another approach towards factorization 
with the help of the Gauss sum, fl56l) . For this technique we assume that we have 
sufficient control over the rescaled chirp ^ to tune it precisely to an integer ^ = I. As 
a consequence, the term linear in the summation index in the phase factor drops out 
and the probability amplitude c^f^\ approximately given by (1571) . is proportional to the 
Gauss sum 



Hence, we deal with a Gauss sum over purely quadratic phases [T5]- 
In Ref. |1] we have analyzed the properties of the function Sn = Sn{(-). In 
particular, we have shown that Sn{^) allows us to factor numbers in a rather 
straightforward way. Since iSiv(£) approximates the excitation probability amphtude 
Ce^^'' of a Floquet ladder given by (H71) . the occupation probability |ce^'^''p at integer 
values I of the dimensionless chirp ^ should yield information about the factors of an 
appropriately encoded number A^. 

In Figure we verify this statement by presenting numerical results for the 
factorization of A^ = 105 = 3 ■ 5 ■ 7 based on (H7I) . In contrast to the previous scheme the 
signal |ci^'^''p is depicted only for integer values of the rescaled chirp ^ = i. Moreover, 
data points with £ being a factor of A^ and their products arrange themselves on a 
straight line through the origin. Data points corresponding to integer multiples of a 
factor are characterized by identical values. On the other hand the signal is suppressed 
at non-factors of A^ in complete accordance with the predictions of Ref. [1] . 




m 



(58) 
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Figure 4. Factorization ofA^ = 21 = 3- 7 with 78 satellites in the Floquet ladder of 
the excited state. Here we show the fluorescence signal |ci''^^^p, (|T7)) . as a function of 
the continuous rescaled chirp ^ for the phase \(p\ = 7r/2 (top) and (p ^ (bottom) of the 
cw-field. The electric field parameters are chosen to yield the width of the weight factor 
distribution An = 12.71 and the modulation index k = 100 • 27r + 7r/4. The positions 
of candidate prime factors are indicated by vertical lines. The insets demonstrate that 
the signal exhibits pronounced maxima at the prime factors £ = 3 and 7 (top) but not 
a.t £ = 2 and 5 (bottom). For the phase (p = the factorization criterion of observing 
distinct maxima at factors of N does not apply. Here the factors are identified by zeros 
(top) rather than maxima. The non- factors have a non- vanishing signal (bottom). 
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Figure 5. Factorization of = 105 = 3 • 5 • 7 in the Floquet-ladder approach using 
the signal Ici^^"*!^, (|T7)) . for integer values ^ — i of the rescaled chirp. At the prime 
factors £ = 3, 5 and 7 and products £ = 15, 21, 35 the signal displays maxima. At 
non-factors the signal is suppressed. Data points corresponding to factors of N are 
situated on a line through the origin. Integer multiples of a factor are characterized by 
the same value of the signal as illustrated by the two horizontal lines.. Since Ce is 
only an approximation of Sm{(-), there are small deviations of this behavior. In order 
to satisfy the criterion of slowly varying weight factors we have chosen the parameters 
An 90 and k = 10^ • 27r + 7r/4. 

6. Pulse train 

In this section we turn to yet another realization of a Gauss sum in a physical system. In 
contrast to the method of the preceding section now the quadratic phase factors are not 
due to a chirped laser pulse, but arise from the combination of a linear time-variation 
of the resonance condition and a pulse train as shown by Figure [61 The probability 
amplitude Ce of excitation after a sequence of laser pulses follows from the sum over 
the contributions from the individual pulses and is of the form of a Gauss sum. Again 
the system is capable of factoring numbers. However, the roles of the trial factor and 
the number to be factored are interchanged. 

6.1. Excitation probability in the weak field limit 

We modulate the energy of the excited state by an electric field 

Em{t)=J^o^, (59) 

which increases linearly in time. Here J-q denotes the amplitude of the field and T is a 
time scale. 

When we substitute this field into the definition (1251) of the phase /3 we find the 
expression 

m = l^t'-/3o^a{t)-/3o (60) 
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Figure 6. Excitation by a pulse train. We consider a two-level system with ground 
state \g) and excited state |e). The modulation field Eui = Em{t), ([59l) . causes a linear 
variation of the excited state energy. Simultaneously a sequence of 2 Af + 1 delta-shaped 
weak laser pulses with envelope h = hit), (|6T1) and carrier frequency cj^ — ujq — 5 
induces a transfer to | e) . Consecutive pulses here depicted for M = 1 are separated by 
T. 



which contains the time-independent Rabi frequency f2ee = pee-^o/^ and (3q = 

^eetl/{2T). 

Moreover, we drive the one-photon transition with the electric field E^, given by 
( 129|) and consisting of a train 

1 ^ 

^^^)^^MTi ^ ^^^-^T) (61) 

~'~ n=-M 

of 2M -|- 1 delta-shaped pulses separated by T. Here, we have chosen a normalization 

oo 

j dt h{t) = 1 (62) 

— oo 

The approximation of the pulse by a delta function reflects the fact that the temporal 
width of the individual pulses has to be small compared to T. 

When we substitute the pulse train h = h(t), (!6T]) . into ( 130|) and perform the 
integration we arrive at 

crHt)=.fi,,e-«^^ Y: exp 

~'~ n=-M 

Here we have assumed that the range of integration in (130|) is large enough to cover 
the whole pulse train. 

Again the probability amplitude Ce of excitation involves the sum over quadratic 
phase factors and is therefore a Gauss sum. 

6.2. Emergence of reciprocate Gauss sum 

Though in principle we could apply the same factorization scheme as in Sec. |5] we 
propose here a more powerful technique for factorization. Indeed, by a proper choice of 



r rr ^eerr 2 

I n 1 n 

2 



(63) 
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parameters we eliminate in f l63p the phase hnear in n. For this purpose we relate the 
detuning 6 and the pulse separation T to the number to be factored by 
6T 

With this choice we find for the quadratic phase 



1 



N 



(64) 
(65) 
(66) 



2 ^ 
where we have introduced the dimensionless variable 
2S_ 

As a consequence, the probability amplitude c^f"^^ for the pulse train given by ([63 
reduces to 

cf^)(t) = ifi,ee'"(*U;v(0 
and is governed by the Gauss sum 

M 



(67) 



n=-M 



-27rin^ 



N 
J 



(6J 



In contrast to the Gauss sums of the preceding sections the roles of A^ and ^ are 
interchanged. Indeed, now the variable ^ appears in the denominator and the number 
A^ to be factored in the numerator. 



6.3. Factorization 

In Sec j5.3l we have shown that the Gauss sum arising in the excitation of the Floquet 
ladder reveals the factors of A^ for a continuous tuning of the chirp ^ as well as for 
integer values ^ = i. Likewise, the Gauss sum An = An{0 defined by ( l68l) provides us 
with the factors for continuous as well as integer values of ^. However, the analysis for 
continuous ^ is more complicated and has been presented in Ref. [16]. For this reason 
we focus in the present section only on the discrete case. 

Since the Rabi frequency Qee is a free parameter we can adjust ^ to be an integer 
£. As a result we arrive at the su 

1 

n=-M 

We now demonstrate that is even more suited to factor numbers than the 
two Gauss sums S or Sn given by fISB]) and (jSHD, respectively. Whenever the integer 
argument £ is a factor q of N the phase of each phase factor of An is an integer multiple 
of 27r. As a consequence, each term in the sum is unity. Since the sum contains 2M + 1 
terms the signal at a factor q of N takes on the maximum value of 

\AN{q)\ = l. (70) 

I In Ref. [Ij we have shown that the Gauss reciprocity relation establishes the connection between the 
two types of Gauss sums and of and respectively. 



-27rin^ 



N 



(69) 
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Figure 7. Factorization of = 1911 = 3 • 7^ • 13 using a train of 21 pulses. The 
modulus |ci^^'| of the signal, (p7)) . exhibits clear maxima indicated by dashed lines at 
integer arguments £ corresponding to factors of N. 



In Figure [7] we illustrate the power of this read-out mechanism of factors using the 
example = 1911 = 3 ■ 7^ ■ 13. We find that already 21 pulses allow us to decide 
whether £ is a factor of N or not. 

7. Comparison of factorization schemes 

We devote this section to a brief comparison of the factorization schemes based on the 
Floquet ladder and the pulse train discussed in SecsO and O Here we first concentrate 
on the methods of readout and then briefly address experimental requirements and the 
necessary resources. 

In the Floquet- ladder approach two techniques to analyze the fluorescence signal 
determined by the population |ce'^^''p in the excited state and given by (H71) . offer 
themselves: (i) We measure |ci^^''p as a function of the continuous chirp ^. In this 
case pronounced maxima at trial factors indicate factors of A^. (ii) An alternative read- 
out relies on the measurement of the signal at integer values = L Here we find that 
the signals at factors of A^ form a straight line through the origin. 

For the pulse-train approach the proposed read-out scheme is based on a 
measurement of the signal |ci^"^^p, (1671) . at integer values of the argument ^. Factors 
of A^ are characterized by the same maximal value, whereas the signal at non-factors is 
suppressed. 

Next we address the experimental requirements for these schemes to work. To 
reveal the factors of a given number A^ it is sufficient to analyze the fluorescence signal 
for values of the dimensionless chirp ^ in the interval [0, a/ZV]. For the continuous version 
the resolution in ^ has to be sufficiently high to resolve the shape of the signal in the 
vicinity of candidate primes. For the discrete scheme the signal has to be acquired only 
for integer arguments I. Nevertheless, we require precise control of ^. When we compare 
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the number of measurements necessary in both schemes to obtain enough information 
for a decision on the factors, the discrete factorization schemes are favorable since less 
data points are required. 

It is also interesting to compare the number of terms in the Gauss sums and 
An necessary to find factors. In the approach based on the Floquet ladder this number 
is determined by the width An of the weight factor distribution Wm, fl55l) . Indeed, this 
distribution has to be sufficiently broad in order to achieve a signal with an appropriate 
contrast. For the pulse-train approach the number of terms contributing to An is 
determined by the number of pulses in the train. Already with a few terms the signal 
has enough contrast to bring out the factors. 

One may wonder whether the elimination of the phase linear in n in the pulse-train 
approach is also possible in the Floquet-ladder system. The basic idea was to chose the 
number N to be factored such that the term linear in the summation index drops out 
of the phase factor. Here, we had three parameters at our disposal: two are required to 
encode the number N and one parameter is free to vary the argument i. 

In the Floquet-ladder approach we also have control over the three parameters 
6, A and 0" for encoding both the number to be factored and the dimensionless 
argument £. However, if all three would have been used to encode N, we would not 
have a parameter left for controlling i. 

8. Conclusions 

In the present article we have proposed three physical systems to implement three types 
of Gauss sums. Our ultimate goal was to construct an analogue computer which would 
calculate these Gauss sums. We have then analyzed the signal to deduce from it the 
factors of an appropriately encoded integer N. 

Our first system is based on a two-photon transition in a ladder system driven by a 
chirped laser pulse. Though this factorization scheme performs well for small numbers 
its performance for larger numbers is questionable since the required dimension D of 
the harmonic ladder needs to be of the order of N. Moreover, it is rather difficult to 
find such an equidistant ladder system in nature. 

For this reason we have investigated two other systems. In the approach of the 
Floquet ladder a cw-field modulates the excited state of a two-level atom giving rise to 
a manifold of equidistant sidebands. When driven by a chirped laser pulse the resulting 
excitation probability amplitude is a Gauss sum. The second technique is based on a 
linear chirp of the excited state energy. A pulse train of delta-shaped pulses ensures 
that the excitation probability amplitude is of the form of a Gauss sum. The origin of 
the quadratic phase factors is different in these two realizations of Gauss sums. In the 
first one they are due to the chirped laser pulse, whereas in the ladder approach they 
originate from a linear chirp of the resonance condition. 

In all three examples the excited state probability is experimentally accessible via a 
detection of the fluorescence signal. Moreover, for each system we have developed rules 



Factorization with Gauss sums: Implementations 



21 



for determining the factors of an appropriately encoded number. 

Our factorization scheme rests solely on interference which implies that the required 
resources grow exponentially with the number of digits of N. This feature is in contrast 
to Shor's algorithm which achieves an exponential speed-up due to entanglement. The 
next challenge is to combine these ideas with entanglement and create a Shor-algorithm 
with Gauss sums. However, this task goes beyond the scope of the present article and 
has to await a future publication. 
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